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Abstract 

In this talk, I present the status of attempts to analyze the behavior of the so-called spatial 
't Hooft loop, which can be taken as an order parameter for the deconfinement phase transition 
in pure SU(iV) gauge theory. While lattice data show a strikingly universal scaling of extracted 
fc-string tensions for various values of k and N, the analytic approach to these observables might 
need some refinement. 
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1. Introduction 

Four-dimensional SU(iV) gauge theory undergoes a phase transition from a cold and 
confined phase into a hot and deconfined phase at a critical temperature T c , which is 
intimately related to the breaking of the Z(N) centergroup symmetry. In this respect, 
a widely used order parameter is the Polyakov loop P = TrPexp(i<7 J^ T dx A ), 
whose thermal average (P) ^ in the deconfined phase. There are however some formal 
issues associated with the use of (P) as order parameter: P cannot be defined at strictly 
zero temperature, and its lattice-discretized version exhibits ultraviolet divergences in 
the continuum limit. 

Analyzing other choices, and noting that the spatial Wilson loop shows only area 
law behavior at T ^ 0, while at T = Q it exhibits area/perimeter law behavior for 
broken/unbroken Z(7V) symmetry, it has been proposed to use the 't Hooft flux loop 
operator [T] - which in some sense is dual to Wilson loops - as an alternative order 
parameter [2]- 

In order to give a definition of the 't Hooft flux loop operator, consider measuring 
color-electric flux through a large surface S, in SU(iV) gauge theory at finite temperature. 
Since the color-electric field (E)j = Poi is in the Lie Algebra of SU(iV), one asks for its 
flux projected onto some special direction Y in SU(iV): $ = - J s dS ■ TrEF. Note that 
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Fig. 1. Comparison of lattice data and analytical results for ratios of the fc-string tensions of Eq. {[2} 
as a function of temperature, for different gauge groups. The horizontal lines are analytic leading-order 
(LO) results the curved lines are NNLO with T c /Aj^ = (1.1, 1.35). The plots are from ref. 0]. 



(<&) = 0. Now define the 't Hooft loop (or E-flux) operator [112] V = exp(47ri$), whose 
thermal average (V) behaves non-trivially. 

The directions (in matrix space) Y are constrained by gauge invariance. Since F^ v = 
VtF^Q, transforms as an adjoint, one has to choose special Y's, namely generators of 
the center of SU(iV) which respect exp(2iriYk) = exp(27ri-^)ljv X Ar for k = 0, ...,N — 1. 
These are not unique; one simple choice [5] is Y k = dia,g({k} N ~ k , {k — N} k ). One hence 
obtains a set of gauge invariant (for a detailed proo{3 see [7]) operators 

V k = exp I — / dS ■ TrEr fe I . (1) 



s 



One can finally define the thermal averages 



(V k ) ~exp(-a k (T)A s ) , (2) 

which at large areas As of the surface S define the fc-string tensions a k . 

The tensions a k can be computed either with analytic methods in a weak-coupling 
expansion (see Sec. ©), or on the lattice 4 J. From Fig. ([T]), one observes fairly good 
agreement between analytic and numerical results for the dimensionless ratios o k ja\. 

Interestingly, when choosing a different normalization on the vertical axis, on the lattice 
one clearly observes more violation of "Casimir scaling" at low temperature [4] , as shown 
in Fig. (2j). Moreover, and most strikingly, the lattice data for different gauge groups and 



1 To prove that Vf, acting on physical states is gauge invariant, the main idea is to construct 
the physical subspace spanned by gauge invariant combinations of gauge fields (which are Wilson 
loops W = TiVexp(i <f c dl ■ A)), to show V k WV^ = exp(2iriY k ) W and therefore V^W(V^ = 
cxp(2ni^lY k ^l~ 1 ) W = Q exp(27riYj c ) Q~ 1 W = exp(2niY k ) W, such that when acting on a physical state 
constructed from one Wilson loop one obtains V^plphys) = VP W (VP)> Vp\0). 



2 



2.2 



2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 



GKA T/A, 



SU(3) 
SU(4), k=l 
SU(4), k=2 
SU(6), k=1 
SU(6), k=2 
SU(6), k=3 
SU(8), k=1 
SU(8), k=2 
SU(8), k=3 
SU(8), k=4 
MSbap 1 ' 35 



-Br- 



1.5 



2.5 



3.5 



4.5 



T/T. 



Fig. 2. The same data as in Fig. JTJ, with a different normalization. Note the universality for different N, 
k. Solid line: NLO o- fe with 2-loop running coupling, at T c /A-^ = 1.35. The plot is taken from ref. [I]. 

different values of k appear to fall onto a universal curve, hinting towards (large-iV?) 
universality. Clearly, it would be nice to understand this behavior analytically. 

2. Analytic results 



The weak-coupling expansion has been driven to next-to-next-to leading order (NNLO), 



<Jk_ _ 'in 2 
T 2 ~ 3v^V 



k(N - k) 1 - 15.2785.. 



-JV 



(4tt) 2 



ca(N,k) 



g 3 N? 
(4tt) 3 



c a {N, k) 



g 4 N 2 



.(3) 



At leading order, one observes Casimir scaling, while the NLO coefficient [5] scales in the 
same way with N and k, provided the expansion is organized in terms of g 2 N . 

The NNLO coefficient C3 (N, k) has been determined in [5] . It is small (which is the 
reason why, for simplicity, in Fig. ^ only NLO is shown); it contributes (in analogy 
to e.g. weak-coupling results for the QCD pressure 0) with a sign opposite to the g 2 
term; its nonanalytic (in a s ) magnitude g 3 is well-understood, originating from infrared 
sensitive 0(g 4 ) graphs that necessitate an all-order resummation of this specific sector; 
and it is expressed in [6] as a specific two-dimensional integral plus an infinite sum. 

From the analytic side, the question now clearly is whether the unknown NNNLO 
term C4(N,k) can be computed, and whether it will help in explaining the apparent 
discrepancy at low temperatures T < 4T C between lattice and analytic results as shown 
in Fig. ([2]). To set our minds for doing this, it might be useful to recall the setup in which 
the analytic results that are known so far have been obtained. 

The basic idea is to compute as tunneling effect through a perturbatively calculable 
potential barrier. Hence, one first needs to compute an effective potential. The presence 
of Vfe breaks the center group symmetry, since it chooses specific directions Yfe. This 
leads to the Polyakov loop P picking up a phase exp(27ri-^) when moving through the 
surface S (at, say, z = 0, if the S is simply chosen as spanned by an L x x L y loop in 
the xy plane), which means that the field component Aq(z) is discontinuous at S. One 



3 



can parameterize this behavior with (diagonal, traceless, N x N) matrices C, according 
to P = Tr exp(iC(z)). Now, constraining the path integral to these allowed values of Aq 
with a delta function one obtains 

(V k )= J[VAc][VA]6(P-Trexp(iC(z)))exp(-S[A]) ~ exp(-L x L y a k ) , 

where the last step is valid in the limit of large loops and defines the minimum value of 
the effective potential U , with a k = mm/c(z)} U[C{z)\. This minimization is taken with 
respect to all C(z) which have a 2nk discontinuity at z — and vanish as z — » oo. 

The minimal profile of the effective potential U is conjectured to be realized along 
the simplest path in the space of the matrices C ~ 2TrqY k , with < q < 1 (while 
there is no general proof yet, this has been shown at large N and for N = 3,4). To 
compute U in a perturbative expansion around the Polyakov loop constraint, one sets 
A = C(z) + gQ and Ai = gQi, chooses background field gauge 5 g f = TTr[£) M (C)(5/i] 2 , 
performs the loop expansion, and finally minimizes with respect to q. On the technical 
side, the g-dependence represents a smooth interpolation between bosonic and fermionic 
Matsubara sums [9]. Additionally, going beyond NLO entails resumming the IR enhanced 
<? 4 diagrams, which is most transparently done in the framework of dimensional reduction. 
As the magnetic sector does not contribute before g 5 , the relative order g A is computable. 

3. Conclusions 

The spatial 't Hooft loop 14 measures the color-electric flux, which at large temper- 
atures is due to free screened gluons. This induces an area law behavior, whose corre- 
sponding tensions a k ~ hi(Vfe) are hard to measure on the lattice, since they are expo- 
nentially suppressed, and hence trigger a sampling problem that increases with the loop 
size. Existing lattice measurements indicate universal behavior of the a k , while NNLO 
weak-coupling results seem to correctly predict ratios of a k , but not their normalization 
close to T c . Hope is expressed that knowledge of the NNNLO term would cure this lack 
of predictivity, and the setup in which to obtain this contribution is roughly sketched. 
As a final remark, let us stress that the ratios o k jo\ can serve as stringent tests for 
formulations of Yang-Mills theory derived from e.g. string theory. 



References 

[1] G. 't Hooft, Nucl. Phys. B 138 (1978) 1. 

[2] A. Kovner and B. Rosenstein, Int. J. Mod. Phys. A 7 (1992) 7419. 
[3] C. Korthals Altes and A. Kovner, Phys. Rev. D 62 (2000) 096008. 
[4] P. de Forcrand, B. Lucini and D. Noth, PoS LAT2005 (2006) 323. 

[5] T. Bhattacharya, A. Gocksch, C. Korthals Altes and R. D. Pisarski, Phys. Rev. Lett. 66 (1991) 998; 

Nucl. Phys. B 383 (1992) 497. 
[6] P. Giovannangeli and C. Korthals Altes, Nucl. Phys. B 721 (2005) 1; Nucl. Phys. B 721 (2005) 25. 
[7] C. Korthals Altes, lecture notes, Bielefeld Graduate School, Sept 2007; see 

http : //www.physik.uni-bielef eld . de/igs/schools/Fall2007/altes_writeup .pdf . 
[8] J. I. Kapust a, Nucl. Phys. B 1 48 (1979) 461. 
[9] Y. Schroder. larXTv : 0807. 0500 [hep-ph]. 



4 



